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Abstract 

It is well known that confining a liquid into a pore strongly alters the liquid behavior. 
Investigations of the effect of confinement are of great importance for many scientific 
and technological applications. Here we present a molecular dynamics study of the 
behavior of cyclohexane confined in carbon slit pores. The local structure and orienta¬ 
tional ordering of cyclohexane molecules are investigated. It is shown that the system 
freezes with decreasing the pore width, and the freezing temperature of nanoconfined 
cyclohenae is higher than the bulk one. 
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I INTRODUCTION 


It is well known that the nature of spatial ordering of a molecular system depends on the 
dimensionality of the space to which it is conhned. Conhnement can modify the dynamics, 
the thermodynamics and the structural properties of liquids when the range of molecular 
interaction and the length scale associated with position correlations in a system are simi¬ 
lar to the length scale of conhnement-. As a result, huids conhned in a pore demonstrate 
the behavior strongly diherent from the bulk case. In slit pores the melting temperature 
depends strongly on the pore width H and the nature of the huid-wall and huid-huid attrac¬ 
tive interaction-. The parameter H determines not only the magnitude of the shift of the 
melting temperature but also the structure of conhned huids and appearance of the possible 
new phases. The density prohle becomes modulated unlike the isotropic one in bulk liquids. 
The study of diherent nanoconhned systems becomes important for both fundamental re¬ 
search and technology because there are a lot of real physical and biological phenomena and 
processes that depend on the properties of such systems. Nano conhnement is considerably 
interesting also due to the new physics observed in these systems. 

In spite of great importance of this issue, there is a lack of systematic investigations of 
ehect of conhnement on the behavior of aromatic hydrocarbons. Only the case of benzene 
received a lot of attention in the literatures^— . However, liquids with more complex chemical 
structure are also of interest. 

In this paper, we present a systematic molecular dynamics study of cyclohexane in the slit 
carbon pores. In our previous publications we considered the behavior of benzene in slit pores 
with graphite and amorphous carbon wallsSi and inside an armchair nanotubeSI. Benzene 
has chemical formula CqHq and planar molecule. The chemical formula of cyclohexane is 
CqH 12- Adding six more hydrogens changes the bonding inside the carbon ring and therefore 
alters the interaction of the ring with external molecules. Moreover, the molecule changes 
its shape: it is not a planar, but warped. In this respect it is interesting to compare the 
behavior of benzene and cyclohexane in the same kinds of pores. Three types of differences 
can affect the results: the different bonding inside the molecules, the different geometric 
shape of the molecules and the different number of hydrogens. 
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The planar structure of the benzene molecule corresponds to the in-plane sp^ hybrid or¬ 
bitals of C, which produce a resonant bond shared by all C atoms. The cyclohexane molecule 
has six sp^ carbon atoms. The strong vr-electron interaction of benzene with nanopore sur¬ 
faces was found-”-*^ which plays an essential role in the formation of the local structure 
of nanoconhned benzene^ii^. At the same time the stronger cyclohexane absorption on 
graphite than benzene absorption was observed—, and one can suppose that the interaction 
of sp^ carbon atoms of cyclohexane with the graphite confining walls can lead to different 
structure and thermodynamic behavior of cyclohexane in comparison with benzene. 

The behavior of cyclohexane in nanopores is interesting in its own rights. In particular, 
in Ref.— it was proposed to use cyclohexane for calibration of NMR experiments which can 
be important for development of experimental techniques. Several publications described 
the behavior of cyclohexane in mica pores^i^^. For example, in Refs.— the crystallization 
of cyclohexane conhned between parallel mica surfaces was studied. It was found that with 
decreasing the pore width up to about Anm (6 or 7 layers of cyclohexane) the crystallization 
occurs which is observed from the sudden increase of a yield stress. It is interesting that it 
was found the remarkable increase of the melting temperature in comparison with the bulk 
case. In experiments^^— the solidihcation was observed at pressures close to the atmospheric 
ones, which gives the possibility to suppose that the geometric confinement alone can lead to 
the solidlike behavior. However, in Ref.— this result was opposed, because in this work the 
experiment was repeated, and no evidence of the increase of melting temperature was found. 
The qualitative support for the results of Refs.— was obtained from computer simulations 
of model systems^i^^. 

An analysis of Ref.— suggests that the sign of the change of the melting temperature 
under confinement depends on the ratio of the fluid-wall to the fluid-fluid attractive inter¬ 
action. A number of studies of different systems confined within activated carbon hbers 
show that the freezing temperature can increase or decrease depending on the nature of the 
system. In Ref.— it was found an increase about 65iF in the freezing temperature of benzene 
in activated carbons. On the other hand, for water, a considerable decrease of the freezing 
temperature was observed— >2^. 

It is interesting to note that in some cases the intermediate hexatic phase can exist— >2^. 
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Figure 1: (a) Cyclohexane molecule. Larger balls represent carbon atoms while the smaller 
ones - hydrogens, (b) A cartoon of the simulation box (see the text). 

As it was proposed in Ref.— the hexatic phase can be observed for the small number of 
layers. For example, in CCI 4 the hexatic phase was found for two layers^. In the case of 
water the hexatic phase does exist in the case of strictly two-dimensional model— which 
demonstrates waterlike behavior—^—. 

In the present paper, we study the local structure and freezing transition of system of 
cyclohexane molecules conhned in graphite slit pores using the realistic AIREBO potential—. 
The molecular models and simulation methods are described in Sec. II, and the results are 
presented in Sec. Ill and discussed in Sec. IV. 


II METHODOLOGY 

We start with brief description of cyclohexane molecule. It consists of a six carbon ring with 
two hydrogen atoms bonded to each of the ring sites (Fig. [1]). Unlike benzene, the carbon 
ring of cyclohexane is not flat, but warped. The total size of the molecule is approximately 

s.sA. 

The simulation methodology is similar to the one used in our previous publication for 
benzeneii. We studied 200 molecules of cyclohexane in slit pores. The walls of the pores 
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consisted of two graphite layers. The inner layer was mobile while the onter one was held 
rigid in order to stabilize the system. The distance between the onter graphite layers is 
denoted as H. The inner layers are separated from the outer ones by 3.4A. We denote 
the distance between the inner layers as h. Obviously, h = [H — 6.8)A. The parameter h 
determines the area in 2 ; direction available for cyclohexane molecules. 

The system was simulated by molecular dynamics method in NVT ensemble (constant 
number of particles, volume and temperature). The temperature was maintained at T = 
SOOit' by means of Nose-Hoover thermostat. The timestep was dtgr = 0.001 fs. The system 
was equilibrated for 1.5 • 10® steps. Then 5.0 ■ 10® steps were done in order to calculate the 
averages. 

In order to see the effect of the walls separation on the structure of the liquid we simulate 
the system at different pore heights. The system sizes in x and y direction are = 49.1902A 
and Ly = 46.86A. The size in z direction (from the bottom of the lower wall up to the top of 
the upper wall) changes from = 19.355A up to if™®* = 33.355A with step dH = 2.0A. 
In total, 8 different pore sizes were studied. 

The system was periodic in XY plane, but confined along axis which is perpendicular 
to the walls. 

All interactions in the system were modelled by AIREBO potential—. This potential 
was specially developed for simulation of hydrocarbon systems. It includes carbon-carbon, 
carbon-hydrogen and hydrogen-hydrogen interactions, i.e. all interactions involved in the 
systems under investigation. 

All simulations reported in this article were done in lammps simulation package^. 

Ill RESULTS 

As many other liquids, cyclohexane in a slit pore forms layered structures. The number 
densities of carbon, hydrogen and centers of mass of the molecules in the case of different 
pore heights are shown in Fig. [2l One can see that in the case of small pores {H < 27.355A) 
three layers are formed: two next to the lower and upper walls and one in the center, while 
at larger separation the peak in the center splits into two ones. It is clearly seen from the 
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density distribution of carbons and centers of mass. At the same time the density of hydrogen 
atoms looks almost flat. 

Below we call the layers close to the upper and lower walls as the outer ones while the 
layer or two layers in the center of the pore as the inner ones. 

The density distribution of the carbon atoms is very sensitive to the pore height. One can 
see that at the height as small as 21.355A a shoulder on the outer peaks appears. Later on it 
transforms into a small subpeak. At H = 27.355A the inner peak splits into two subpeaks. 
On increasing H these subpeaks become separate peaks. 

The density distribution of CM are less sensitive to the change of H. One can not see 
small subpeaks and splitting aX H = 27.355A. However, one observes that the central peak 
splits into two well dehned peaks on change from H = 29.355A to H = 33.355A. 

Next we consider the orientation of the cyclohexane molecules inside the pore. It is 
common to characterize the orientation of molecules by second order Legender polynomial 
P 2 {cos{ 6 )) = 1.5cos^{6) —0.5, where 6 is the angle between some vector which can determine 
the molecule orientation (for example, in the case of benzene it is the normal vector to the 
plane of the benzene ring) and z axis. P 2 {cos{ 6 )) = 1 means that the molecule is parallel to 
the wall, while P 2 {cos{ 6 )) = —0.5 corresponds to perpendicular orientation of the molecule. 
In our previous publication we considered benzene confined in graphite por©^. It was shown 
that in the case of small pores the molecules lie parallel to the walls forming a planar 
structure. However, it can be related to the planar shape of the benzene molecules. Unlike 
the case of benzene the molecules of cyclohexane are not planar, which makes the situation 
more complex. First of all, in the case of benzene one characterize the orientation of the 
molecule by the normal vector to the benzene ring plane. The carbon ring of cyclohexane 
molecule is not planar so one needs to adjust the dehnition. 

A procedure for non-planar molecules was proposed in our previous publication devoted 
to the benzene inside a carbon nanotub©^. For the sake of completeness we give a brief 
description of this procedure below. 

Let us denote all carbon atoms in the ring by numbers from 1 to 6 and define the vectors 
connecting these atoms; ri 2 , ris and so on. Each pair of such vectors can be used to dehne 
a plane. The normal vector to this plane can be calculated as the vector product of the 
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Figure 2: Distribution of number density of the molecules centers of mass (CM), carbon (C) 

and hydrogen (H) atoms perpendicular to the graphite wall for several pore sizes. In order 

to have the same scale we multiply the CM curves by 6 and divide the H curves by 2. The 

sizes of the pores are shown in the corresponding plots. 
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Figure 3: The distribution of second order Legender polynomial P 2 {cos{ 6 )) along z axis of 
the system for different pore sizes. 





































































corresponding vectors between the atoms. Employing this procedure we can calculate 20 
vectors perpendicular to the planes of triplets of the carbon atoms and calculate P2{cos{9)), 
where 6 is the angle between the 2 ; axis and the normal vector of each plane. The final 
orientational parameter P 2 is dehned as an arithmetic average of all 20 values for different 
triplets. 

The distributions of P 2 along the z axis are shown in Figs. |3] (a)-(g). One can see that 
in all considered cases the molecules at the outer layers lye almost parallel to the graphite 
sheets. At PI = 19.355A the corresponding peaks are almost unity. Upon increasing of H 
these peaks slightly decay, but even at H = 33.355A they are still as high as 0.7. 

Unlike the outer peaks the inner one rapidly smashes out. One can see that already at 
H = 25 . 355 A it practically disappears. Later on it splits into two peaks, but both of them 
are small and noisy which means that the orientational order is almost lost. 

Importantly, in between of the layers the most preferable orientation of the molecules is 
mainly perpendicular to the wall. It means that the molecules located between the layers 
construct some kind of bridges between them. This phenomena can be important for stability 
of the layered structure. 

As it was described above the system can be divided into layers. These layers can be 
considered as quasi two dimensional systems. This allows us to consider in-plane structure 
of each layer. The structure of these two-dimensional layers can be characterized by in-plane 
radial distribution functions (72 (^)- 

In order to find the borders of the layers we took the minima of the density distribution 
of the centers of mass (Figs. [2]). These borders are given in Table II. Note that the width 
of the layers can be as large as 5 — 6 A, i.e. the layers are not perfectly flat. Therefore 
the two-dimensional representation is not exact. However, since the width of the layer does 
not exceed the size of the molecules quasi two-dimensional consideration can be a useful 
approximate tool. 

Two-dimensional radial distribution functions of centers of mass of the molecules in dif¬ 
ferent layers are shown in Figs. |l](a)-(h). The results for outer peaks of p{z) (Fig. ??)(curve 
’wall’) and for inner ones (curve ’center’) are shown. For the outer peaks of p{z) the radial 
distribution function g 2 {r) is averaged over the upper and the lower layers. The same aver- 
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H=19.355 A 








Figure 4: In-plane radial distribution functions of centers of mass of CqHu molecules in 
graphite pores. The data are averaged for two peaks next to the wall (’wall’ curve in the 
plots). In the case of the peaks at the center at H > 29.355A the curve is averaged over two 
central peaks. 
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Figure 5: Left panel; Snapshot of the outer layer of centers of mass of CqHi 2 molecules in 
H = 25 . 355 A slit pore. Right panel: the same for the inner layer. 
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Figure 6: Left panel; Snapshot of the outer layer of centers of mass of C^Hu molecules in 
H = 33 . 355 A slit pore. Right panel: the same for the inner layer. 
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aging is applied to the inner layers if the central peak splits into two ones. One can see that 
aX H = I 9 . 355 A the first peak of g 2 {'r) for the central layer is much larger then for the outer 
ones. We suppose that at such low walls separation the dynamics of liquid is too slow and 
the results strongly depend on the initial conhguration. The considerable splitting of the 
second peak of the radial distribution function allows to suppose that in this case there is a 
glassy state of the system, however, more thorough investigation of this point is necessary. 
For H = 25 . 355 A one can observe 4 peaks of g 2 which corresponds to very structured liquid. 
Moreover the positions of the peaks is (for H = 25.355A) ri = 6.23A, r 2 = 11.lA, = 16.9A 
and r 4 = 22.26A. The ratio of r 2 /ri and r^/ri corresponds to hexagonal two-dimensional 
lattice. Fig. [S] shows a snapshot of outer layer (left panel) and inner layer (right panel) for 
this value of H. One can see that in both cases defected triangular lattices are formed. 
One can conclude from both ( 72 (^") for H = 25.355A (Fig. H]) and the snapshots of the layers 
(Fig. [5]) that the inner and outer layers of this system are structurally equivalent. 

The same conclusion holds true for the system with H = 27.355A. In the case of if = 
29 . 355 A, H = 3 I. 355 A, and H = 33.355A the situation looks different. The crystallinity 
of both inner and outer layers drops. This effect is especially noticeable in the inner layers 
where only first peak of g 2 is well pronounced, i.e. the centers of mass of the molecules in 
the inner layers behave as moderate density liquid. g 2 {'f) of outer layer looks like in the case 
of high density liquid. Snapshots of these cases are shown in Fig. [6l 

To make the picture more clear, let us consider the behavior of the transverse and 
longitudinal D^y diffusion coefficients of the system as functions of H. and D^y are 
calculated with the help of the Einstein equation—. From Fig. [3 one can see, that for 
H < 27 . 355 A diffusion coefficient becomes equal to zero. Taking into account Figs. [5] and 
El this behavior can be interpreted as the freezing transition which occurs with decreasing 
H. It should be noted that at 77 = 27.355A the number of layers changes from 4 to 3 
(Fig. [2]). It is in qualitative agreement with the experimental results^^— on the behavior of 
cyclohexane in mica slit pores. Moreover, in our case the freezing temperature is T = 30077 
whereas the bulk freezing temperature of cyclohexane is equal to Tj, = 279.6277. This 
increase of freezing temperature is again in qualitative agreement with the experiment—^— 
and computer simulations of model system&2^>2^. 
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Figure 7: Transverse and longitudinal D^y diffusion coefficients of the system as fnnctions 
of 77 at T = 30077. 


77, A 

h, A 

p, g/cm? 

17.172 

10.462 

1.075 

18.172 

11.462 

0.982 

19.172 

12.462 

0.903 

20.172 

13.462 

0.836 

21.172 

14.462 

0.778 

22.172 

15.462 

0.728 


Table 1: The size of the graphite pore 77, the distance between the inner graphite sheets h 
and the density of conhned benzene p 
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H 

Layer 1 

Layer 2 

Layer 3 

H = 17.172 

= 5.29 

= 7.97 

= 10.68 


Z2 = 6.49 

^2 = 9.17 

^2 = 11.71 

H = 18.172 

= 5.34 

= 8.21 

= 11.19 


Z2 = 6.61 

^2 = 9.7 

Z2 = 12.68 

H = 19.172 

= 5.71 

= 8.78 

= 11.69 


Z2 = 7.01 

Z2 = 10.47 

Z2 = 13.61 

H = 20.172 

zi = 5.77 

zi = 8.96 

zi = 12.43 


Z2 = 7.63 

Z2 = 11.09 

to 

oo 

H = 21.172 

zi = 5.25 

Zi = 9.14 

Zi = 12.03 


Z2 = 9.14 

Z2 = 12.03 

Z2 = 15.16 

H = 22.172 

zi = 4.74 

zi = 9.13 

zi = 12.90 


Z2 = 9.13 

^2 = 12.90 

^2 = 16.81 


Table 2: The coordinates of the borders of benzene layers for different graphite pore size. 

IV CONCLUSIONS 

In conclusion, we present the molecular simulation study in NVT ensemble of the cyclohexane 
confined between the carbon plates (see Fig. [1]) using the realistic AIREBO potential. The 
local structure of the system is studied as a function of the distance between the plates H and 
it is shown that the molecules form the layered structure. The number of layers increases with 
increasing the width of the pore. The orientation of the cyclohexane molecules is also studied. 
It is found that the molecules at the outer layers lye almost parallel to the graphite sheets, 
while in between of the layers the most preferable orientation of the molecules is mainly 
perpendicular to the wall. This behavior is similar to the one of the benzene molecules in 
carbon nanotubes^. 

In our previous publication— we considered benzene molecules in the same kind of 
graphite slit pores. It was shown that at small H the outer benzene layers crystallize in 
a planar structure while the inner one stays relatively disordered. Unlike this, all layers of 
cyclohexane behave qualitatively similar: all layers simultaneously demonstrate freezing or 


14 












melting transitions. However, at large pore sizes H the outer layers become more ordered 
what is intuitively clear: the presence of a crystalline wall with the same symmetry as the 
carbon ring of cyclohexane acts as a periodic external potential on the liquid and makes it 
more ordered. 

From the local structure (Figs. |5]and|6]) and behavior of the diffusion coefficient (Fig. 
one can conclude that the system freezes with decreasing the width of the pore and cor¬ 
responding change of the number of layers from 4 to 3. It is interesting that the freezing 
temperature of the confined cyclohexane is higher than the bulk freezing temperature in qual¬ 
itative agreement with the experiment—^— on cyclohexane conhned in the mica pores. It 
should be noted that no signs of the hexatic phase were found in the present study, however, 
further investigations in this direction are necessary. 

ACKNOWLEDGMENTS 

The authors are grateful to S.M. Stishov and V.V. Brazhkin for valuable discussions. Sim¬ 
ulations were made at supercomputing resources in NRC ’’Kurchatov Institute”, which are 
supported as the center for collective usage (project RFMEFI62114X0006, funded by Min¬ 
istry of Science and Education of Russia). The work was supported by the Russian Science 
Foundation (Grant No 14-12-00820). 


15 


References 


1. G. A. Mansoori, S. A. Rice, Adv. Chem. Phys. 156, 197 (2015). 

2. R. Radhakrishnan, K.E. Gubbins, M. Sliwinska-Bartkowiak, J. chem. Phys. 116 , 1147 

( 2002 ). 

3. L.Z. Zhang, P. Gheng, Phys. Ghem. Gomm. 6, 62 (2003). 

4. M.D. Ventra, S.T. Pantelides, N.D. Lang, Appl. Phys. Lett. 76, 3448 (2000). 

5. Y.G. Ghoi, W.Y. Kim, K.-S. Park, P. Tarakeshwar, K.S. Kim, J. Ghem. Phys. 122 , 
094706 (2005). 

6. E. Joselevich, GhemPhysGhem 5, 619 (2004). 

7. L. M. Woods, S. G. Badescu and T. L. Reinecke, Phys. Rev. B 75, 155415 (2007). 

8. G.Wei, D. Lin and Z. Dongqiang, Environ. Sci. Technol. 41, 8295 (2007). 

9. Th. Tran-Duc, N. Thamwattana and J. M. Hill, J. Math. Ghem. 49, 1115 (2011). 

10. Y. Shim, Y.-J. Jung, and H. J. Kim, Phys. Ghem. Ghem. Phys. 13, 3969 (2011). 

11. Yu. D. Fomin, J. Gomp. Ghem., 34, 2615 (2013). 

12. Yu. D. Fomin, E. N. Tsiok, and V. N. Ryzhov, J. Gomp. Ghem. 36 , 901 (2015), DOL 
10.1002/jcc.23872. 

13. G. Dosseh, Y. Xia, and G. Alba-Simionesco, J. Phys. Ghem. B 107 , 6445 (2003). 

14. R. N. Smith, G. Pierce, and H. Gordes, J. of American Ghemical Society 72 , 5595 (1950). 

15. D. W. Aksnes and L. Kimtys, Appl. Magn. Reson. 23, 51 (2002). 

16. J E. Gurry, H. K. Ghristenson, Langmuir 12, 5729 (1996). 

17. M. Heuberger and M. Zach, Langmuir 19 , 1943 (2003). 

18. J. Klein and E. Kumacheva, Science 269 , 816 (1995). 


16 


19. J. Klein and Eu. Kuniacheva, J. Cheni. Phys. 108 , 6996 (1998). 

20. E. Kuniacheva and J. Klein, J. Chem. Phys. 108, 7010 (1998). 

21. H. Christenson, Colloids Surf. A 123 , 355 (1997). 

22. S. T. Cui, P. T. Cummings, and H. D. Cochran, J. Chem. Phys. 114 , 7189 (2001). 

23. R. Radhakrishnan, K. E. Gubbins, and M. Sliwinska-Bartkowiak, J. Chem. Phys. 116 , 
1147 (2002). 

24. A. Watanabe and K. Kaneko, Chem. Phys. Lett. 305 , 71 (1999). 

25. M. Sliwinska-Bartkowiak, G. Dudziak, R. Sikorski, R. Gras, K. E. Gubbins, and R. 
Radhakrishnan, Phys. Chem. Chem. Phys. 3, 1179 (2001). 

26. S. Han, M. Y. Choi, P. Kumar and H. E. Stanley, Nature Physics 6, 685 (2010). 

27. D.E. Dudalov, Yu.D. Fomin, E.N. Tsiok, and V.N. Ryzhov, J. Phys.: Conference Series 
510 , 012016 (2014) (doi:10.1088/1742-6596/510/l/012016). 

28. D.E. Dudalov, Yu.D. Fomin, E.N. Tsiok, and V.N. Ryzhov, J. Chem. Phys. bf 141, 
18C522 (2014). 

29. D.E. Dudalov, Yu.D. Fomin, E.N. Tsiok, and V.N. Ryzhov, Soft Matter 10 , 4966 (2014). 

30. Y.D. Fomin, N.V. Gribova, V.N. Ryzhov, S.M Stishov, and D. Frenkel, J. Chem. Phys. 
129 , 064512 (2008). 

31. N.V. Gribova, Y.D. Fomin, D. Frenkel, and V.N. Ryzhov, Phys. Rev. E 79 051202 
(2009). 

32. Y.D. Fomin, E.N. Tsiok, and V.N. Ryzhov, J. Chem. Phys. 135 , 234502 (2011). 

33. St.J. Stuart, A.B. Tutein and J.A. Harrison, 112 , 6472 (2000). 

34. S. Plimpton, Fast Parallel Algorithms for Short-Range Molecular Dynamics, J. Comp. 
Phys., 117 , 1 (1995); http://lammps.sandia.gov/ 


17 



35. Daan Frenkel and Berend Smit, Understanding molecular simulation (From Algorithms 
to Applications), 2nd Edition (Academic Press), 2002. 


18 


